Electron Beam Supercollimation in Graphene Superlattices 
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Although electrons and photons are intrinsically different, importing useful concepts in optics to 
electronics performing similar functions has been actively pursued over the last two decades. In par- 
ticular, coUimation of an electron beam is a long-standing goal. We show that ballistic propagation 
of an electron beam with virtual no spatial spreading or diffraction, without a waveguide or external 
magnetic field, can be achieved in graphene under an appropriate class of experimentally feasible 
one-dimensional external periodic potentials. The novel chiral quasi-one-dimensional metallic state 
that the charge carriers are in originates from a collapse of the intrinsic helical nature of the charge 
carriers in graphene owing to the superlattice potential. Beyond providing a new way to construct- 
ing chiral one-dimensional states in two dimensions, our findings should be useful in graphene-based 
electronic devices (e. g. , for information processing) utilizing some of the highly developed concepts 
in optics. 
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Electronic analogues of many optical behaviors such as 
focusing ^, 2r, 3] J coUimation [^] , and interference [5] have 
been achieved in two-dimensional electron gas (2DEG), 
enabling the system as a basic platform to study funda- 
mental problems in quantum mechanics 0, 0, [1] as well as 
quantum information processing The close relation- 
ship between optics and electronics is been made possi- 
ble due to the ballistic transport properties of a high- 
mobility 2DEG created in semiconductor heterostruc- 
tures [10']. Among those electronics-optics analogues, the 
coUimation or quasi-one-dimensional (quasi-lD) motion 
of electrons and photons are particularly important not 
only to achieve electronic quantum devices [7j, |9| but also 
to realize ultracompact integrated light circuits 11, 1^ . 



Usually, electrons originated from a point source may 
be controlled by electrostatics or geometrical constric- 
tions [^,i3i]- Quasi-ID electronic states and focusing have 
been achieved in a 2DEG with the help of external mag- 
netic fields, e. g. , employing magnetic focusing [l[ and 
quantum Hall edge states [3, [l3| ■ However, it would be 
difficult to integrate them into a single electronic device 
due to the external high magnetic field apparatus needed. 
In view of recent successful demonstrations of extreme 
anisotropic light propagation without diffraction, called 
supercollimation in photonic crystals [111, lHj Oj Ell i an 
analogue of this effect in two-dimensional (2D) electron 
systems may also be possible. In this work, we demon- 
strate that graphene 



17l . [18| in an external periodic 



potentials, or a graphene superlattice, is particularly suit- 
able to realize electron supercollimation in two dimen- 
sions. 

The isolation of graphene [H, 17, [I3|, a single layer 
of carbon atoms in a honeycomb structure composed 
of two equivalent sublattices, offers a new dimension to 
study electronics-optics analogues. Carriers in graphene 
exhibit ballistic tra nsp ort on the submicron scale at 
room temperature [iQl and with mobility up to 2 x 



10^cm^V~^s~^ [2^. Graphene electronic states have 
an internal quantum number, a pseudospin, that is not 
found in normal electronic systems and strongly influ- 
ences the dynamics of the charge carriers. The pseu- 
dospin is of central importance to many of the novel phys- 
ical properties of graphene [H, [13, [H, IH, S [Hi Q , and 
it also plays a significant role in the present work. 

The low-energy quasiparticles in graphene whose 
wavevectors are close to the Dirac point K in the Bril- 
louin zone are described by a 2 x 2 Hamiltonian matrix 
i/o(k) = hvQ ( <^xkx ~^ (^yky), where vq w 10*' m/s is the 
band velocity, k is the wavevector measured from the K 
point, and cr's are the Pauli matrices. The energy eigen- 
values are given by (k) — shvo |k| where s = (— 1) 
denotes the conical conduction (valence) band (Fig. IB). 
The sublattice degree of freedom of the quasiparticles 
in graphene can conveniently be described with a pseu- 
dospin basis, or spinors, where the |t) and the |J,) pseu- 
dospin states of (Jz represent 7r-electron orbital on the A 
and B sublattices of the structure of graphene, respec- 
tively. This Hamiltonian is very similar to the one used 
to model neutrinos as massless Dirac fermions 2^ 24 1. 
The corresponding wavefunction is given by 



(1) 



where 9k is the angle of the wavevector k with re- 
spect to the X-axis. Equation ([T]) may be viewed as 
having the pseudospin vector being parallel and anti- 
parallel to the wavevector k for electronic states in the 
upper (s — 1) and the lower (s — —1) band, respec- 
tively (Fig. IB) [P 



24| . As the spin plays a role 



in the dynamics of neutrinos, the present pseudospin 
is similarly irnportant in the quasiparticle dynamics of 
graphene [23], [2J| . 

Now, let us consider a ID external periodic potential 
V{x) applied to graphene. The potential is taken to vary 



2 




FIG. 1: Electron energy dispersion relation in a special graphene superlattice. (A) Schematic diagram of a Kronnig-Penney 
type of potential applied to graphene with strength Uo inside the gray regions and zero outside. The lattice period is L and the 
barrier width is w. (B) Schematic diagram showing the electronic energy dispersion relations and pseudospin vectors (black 
arrows) in graphene. (C) Contour plot of the first electronic band above the Dirac point energy in pristine graphene. The 
energy difference between neighbouring contours is 25 meV, with the lowest contour near the origin having a value of 25 meV. 
(D) The electronic energy dispersion relation E versus kx with fixed ky. Red, green and blue lines correspond to ky = 0, 
0.1 n/L and 0.2 n/L, respectively, as indicated in C. (E), (F) and (G) Same quantities as in B, C and D for the considered 
SGS {Uo ~ 0.72 eV, L = 10 nm and w — 5nm). Red and blue arrows in E represent the 'right' and the 'left' pseudospin state, 
respectively. 



much more slowly than the carbon-carbon distance so 
that inter- valley scattering can be neglected [H,!!!. Un- 
der this condition and for low-energy quasiparticle states 
whose wavevectors are close to the K point, the Hamil- 
tonian reads 



H = hvQ {—iUxOx + o'yky + I V{x)/hvQ) . 



(2) 



where / is the 2x2 identity matrix. The eigenstates and 
eigenenergies of the Hamiltonian H in Eq. ^ may be 
obtained numerically in the general case or analytically 
for small k. 

It has been predicted that, in a graphene superlat- 
tice with a slowly varying ID periodic potential or a 2D 
periodic potential of rectangular symmetry, the group 
velocity of its low-energy charge carriers is renormal- 
ized anisotropically [25| . Unlike bare graphene which 



has an isotropic (zero mass) relativistic energy disper- 
sion (Figs. IB, IC and ID), graphene under some specific 
superlattice potentials displays extremely anisotropic 
quasiparticle energy dispersion: the group velocity near 
the Dirac point along the direction perpendicular to the 
periodicity of the potential vanishes while the one parallel 
to the periodicity direction is intact (25j . 

We consider a Kronig-Penney type of potential with 
barrier height Uq, lattice period L, and barrier width w, 
periodic along the x direction (Fig. lA). These poten- 
tial parameters can be tuned so that the group velocity 
of the quasiparticles (with wavevector close to the Dirac 
point) along the y direction vanishes [2^. We shall focus 
on a graphene superlattice under one of these conditions 
(Uo = 0.72 eV, L = 10 nm, and w = 5 nm) The pa- 
rameters used here are experimentally feasible as shown 
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FIG. 2: Pseudospin collapse in a special graphene super- 
lattice. (A) and (B) Calculated overlap of two quasiparti- 
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cle states i/'3,k(r) and V's'.k'(r), 
graphene versus 6]^ and B^^i which are the angles between the 
X axis and wavevectors k and k' (jkj = |k'| — O.lvr/L), re- 
spectively. The overlap is shown in a gray scale (0 in black 
and 1 in white). The two states are in the same band (s' = s) 
in A and are in different bands (s' = — s) in B. (C) and 
(D) Same quantities as in A and B for the considered SGS 
{Uo = 0.72 eV, L = 10 nm and w = 5nm). 



in recent studies 2^, 27, 2^ 2^. Later, we will relax the 
special condition to confirm the robustness of the pre- 
dicted supercoUimation. 

The quasiparticle energy dispersion of this superlattice 
(Figs. IE, IF and IG) shows that, not only the group ve- 
locity of quasiparticles at the Dirac point along the y 
direction vanishes, there is hardly any dispersion along 
the ky direction within a good fraction of the supercell 
Brillouin zone (Figs. IF and IG). This portion of the en- 
ergy dispersion in this superlattice thus is well described 
by the relation 



(3) 



The deviation of the actual energy dispersion from that 
of Eq. ([3]) is less than 5% for k vector as large as 40% of 
the supercell Brillouin zone considered (Figs. IF and IG). 
This is an equation for wedges. Thus, for some specific 
superlattice potentials, graphene turns from a zero-gap 
semiconductor into a quasi-lD metal with a finite and 
constant density of states about the Dirac point energy. 
We shall call this class of graphene superlattices as special 
graphene superlattices (SGSs). 

Along with the quasiparticle energy dispersion, the in- 
ternal pseudospin symmetry of the electronic states in 



SGSs also undergoes a dramatic alteration. We calculate 

2 



numerically the overlap 



quasiparticle states ips,'k{^) and V's',k'(r). If there were no 
external periodic potential, this overlap would be simply 
(1 4- ss' cos (6*1^' " ^'k))/2 as seen from Eq. ^ (Figs. 2A 
and 2B). The same overlap in the SGS (Figs. 2C and 2D) 
is however dramatically different from that in graphene, 

and can be well described by (ips' ,k'\e''^^ ^^'^'^\ips,kj = 
(1 -I- ss' sgn (kx) sgn {k'^))/2 . This behavior is robust for 
a wide range of the magnitudes of k and k', extending 
over a good fraction of the supercell Brillouin zone with 
a high degree of accuracy. The eigenfunctions of an SGS, 
for states with small k, can be deduced from this result, 
together with results on the numerically obtained wave- 
functions, (also from analytic calculation: see Supporting 
Information) as having the form of 
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^ \^ s sgn (kx) 



(4) 



where f{x) is a real function. Thus, the spinor in Eq. (|4]) 
is an eigenstate of ax- Therefore, the direction of the 
pseudospin is quantized so that it is cither parallel or 
anti-parallel to the x direction, which is the direction of 
the periodicity of the superlattice potential (Fig. IE), and 
not to the wavevector k as is the case in pristine graphene 
(Fig. IB). In other words, the pseudospin in the SGS col- 
lapses into a backward ('left') or a forward ('right') state. 
The resulting quasi-one-dimensionality in the energy dis- 
persion relation and in the pseudospin of quasiparticles 
in the SGS significantly changes the already unique prop- 
erties of graphene. 

The quasi-one-dimensionality and specific chiral na- 
ture of the SGS makes it a natural candidate for elec- 
tron supercoUimation. It is indeed the case that, when a 
wavepacket of electron is injected into an SGS, the prop- 
agating packet exhibits essentially no spatial spreading, 
i.e. , electron beam supercoUimation is realized (Fig. 3). 
We have calculated the time-evolution of a gaussian wave 
packet (with spatial extent along the x and the y direc- 
tion given by 2ax — 40 nm and 2ay — 200 nm, respec- 
tively) composed of states in the first band above the 
Dirac point energy with a central wavevector k^ (Fig. 3). 
To provide a measure of the electron beam coUimation, 
we compute the angle 9c in which direction the beam 
intensity is maximum and the angular spread AO which 
gives half the maximum intensity when the angle is at 
6c ± A0. For a central wavevector kc parallel to the x di- 
rection with energy i?(kc) = Eq = hvo O.Itt/L = 0.02 eV, 
the angular spread A^^ in pristine graphene is A6 = 55 ° 
(Fig. 3A), whereas in the SGS, AO = 0.3 ° (Fig. 3B), 
about 200 times smaller than in graphene. Specifically, 
in the SGS, the spread of the wave packet in the y di- 
rection after proceeding 0.1 mm along the x direction is 
only 500 nm. Therefore, supercoUimation of currents of 
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FIG. 3: Special graphene superlattice as an electron supercoUimator. (A), (E) and (G) Time-integrated probability density 
of an electron wave packet, \'^{x,y,t)\^dt, in graphene. The initial {t = 0) wave packet is a Gaussian localized at the 

coordinates origin (middle of the left edge of each panel) [^^(a;, i/, 0)^ ~ exp [— (^x^/2a'^ + i/^/2(Ty)] where 2crx = 200 nm 
and 2ay — 40 nm. The wave packet in wavevector space is set to be localized around a specific kc. In A, kc is set by 
_B(kc) = £"0 = kvo 0-ln/L = 0.02 eV and kc = x. In E, kc is set by i5(kc) — 2Eo and kc = x. In G, k^ is set by E{\s.c) = Eq 
and kc = x j \f2 -\- ij j \f2. Qc denotes the angle (defined with respect to the a;-direction) along which direction the intensity is 
maximal and AS denotes the angular spread which gives half the maximum intensity when the angle is at Qc ± AS. (B), (F) 
and (H) Same quantities as in A, E and G for the considered SGS (C/o = 0.72 eV, L = 10 nm and w = 5 nm), respectively. 
(C) and (D), Same quantities as in B for graphene superlattices corresponding to a superlattice potential that is otherwise the 
same as the SGS studied but with a period L change of t^LjL = 5% and AL/ L = —5%, respectively. 



nanoscale width in the SGS can, in principle, be achieved 
and maintained as long as the ballistic transport occurs 
in the system. 

Even when the experimental situation deviates from 
the ideal conditions for SGSs, supercoUimation persists. 
Hence the phenomenon is quite robust. First, for ex- 
ample, if we consider an imperfection in making a su- 
perlattice potential such that the periodicity is slightly 
larger or smaller (AL/L = ±5%), the calculated time- 
evolution of a gaussian wave packet shows the angular 
spread = 0.007° and 0.5°, respectively (Figs. 3C and 
3D). Second, when considering doped SGSs or high en- 
ergy electron injection such that E{\s.c) is doubled, the 
angular spread is still very small (Fig. 3F). Third, even 
when kc is off from the coUimation direction {-\-x) by 



45 °, the angular deviation is still negligible {6c = —1.1 °) 
(Fig. 3H). In graphene, on the other hand, the propaga- 
tion direction and spread of the wave packet sensitively 
depends on the magnitude and the direction of the cen- 
tral wavevector kc (Figs. 3E and 3G). This robustness 
here is quite contrary to the case in optics where the ef- 
ficiency of supercoUimator, superprism or superlens (30j 
depends sensitively on the magnitude and the direction 
of the light wavevector and in gen eral pro vides a very 
narrow effective bandwidth ll|, [ij, [l5|. From our 
calculations, we expect that the predicted supercoUima- 
tion be observable in SGSs over a wide operation range. 

Lastly, we consider the tunneling properties of injected 
electrons into SGSs from pristine graphene, which pro- 
vides another measure of the efficiency of electronic de- 
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FIG. 4: Reflection and transmission at a graphene — spe- 
cial graphene superlattice interface and virtual imaging. (A) 
Schematic diagram showing the incident, the reflected and 
the transmitted wave (the band index is set to s = 1) at a 
graphene — SGS interface, with the relative amplitudes being 
1, r and t, respectively. Thick arrows represent the wavevec- 
tors of the corresponding waves. The incidence and reflection 
angle is 6. (B) Reflectance R = |rp versus the incidence 
angle 6. (C) Schematic diagram showing the propagation 
of electron waves in graphene — SGS — graphene geometry. 
Thickness of each arrow is proportional to the actual inten- 
sity of the wave. A virtual image (dashed disk) is formed at 
a place far from the actual wave source (solid disk). 



vices based on SGSs. When an electron is injected into 
an SGS with an incidence angle 9 from the graphene 
side, the wavevector of incident electron is given by 
ki = ko cos 9 x + ko sin 9 y and those of the reflected and 
transmitted electrons by = — fco cos 9 x + ko sin 9 y 
and kt = ko X + ko sin 9 y, respectively (Fig. 4A). Here, 
we have made use of the continuity of the transverse com- 
ponent of the wavevector and conservation of energy, to- 
gether with the novel dispersion relation given by Eq. ^ . 
Using the continuity of the wavefunction in the system 
described by Eqs. H]) and (g]), we find that the reflectance 
R ^ IrP is 



R{9) = tan^ - 



(5) 



Interestingly, the reflectance is independent of the spe- 
cific form of the external periodic potential of the SGS. 
Equation (O indicates that the transmittance is large for 
most incidence angles. For example, even at = 45 °, 
the reflectance is less than 20 % (Fig. 4B). Therefore, 
the SGS is not only an excellent electron supercoUimator 
but also a good transmitter in a graphene-SGS-graphene 
junction (Fig. 4C). Utilizing this property, an immedi- 



ate application could be made to demonstrate an elec- 
tronic analogue of virtual imaging in this configuration 
(Fig. 4C). 

Given the recent rapid progress in graphene superlat- 
tices fabrication 2^ 22, 2^ 2^, the manipulation of elec- 
trons in ways similar to that of photons in optics by using 
the supercoUimation effect discussed here together with 
other optics analogues 2l|, |22| is expected to soon be 
practicable. The SGSs have the promise of playing a 
unique role in devices based on the synergetic importing 
of concepts and techniques well developed in optics to 
electronics. 
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Supporting Online Material: Analytical solution of 
electronic states of special graphene superlattices 



When graphene is in a ID external periodic potential 
V{x), the Hamiltonian for electrons whose wavevectors 
are close to the K pointreads 



H = hvQ {-iuxdx + (Tyky 



I V{x)/hv^) 



(6) 



where / is the 2x2 identity matrix. In this case, we 
assume that the potential varies much more slowly than 
the carbon-carbon distance. So, the inter-valley scatter- 
ing can be neglected safely (24. 31|. 

If a transform H' = HU is applied to Eq. ^ 
with the unitary matrix U — exp (— icTj^ a(x)/2), where 
a{x) = 2 V{x')dx' /hvo (we assume that a constant 
is subtracted from V{x) to set its average to zero), the 
resulting Hamiltonian H' reads 

H' = hvo {~i<Jxdx + ( cosa(a;) ay — sina(x) cTz ) ky) . 

_ (7) 

(A similar transform was applied to the Hamiltonian of 
a carbon nanotube under a sinusoidal potential for the 
specific case of finding the band gap openin g b ehavior 
at the supercell Brillouin zone boundary S^, l^l-) The 
terms having ky could be treated as a perturbation if ky 
is small. The eigenstate of the unperturbed Hamiltonian 
is given by 



^\ssgn {kx) 



(8) 



Within second order perturbation theory, using Eqs. (O 
and ([H]), the energy eigenvalue of the superlattice is 



s',G 



s\kx\ - s'\kx + Gx\ 



(9) 



and the wavefunction is 



s',G 



%s',u+G |cosa(a:) ay - sma{x) a,,\%^ 



s\kJ 



GJ 



k+G 



(r) 



(10) 



Here, G's are the superlattice reciprocal lattice vectors G = mGo, where Go = (27r/i,0) and m is an integer. 
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Therefore, in order for the superlattice to be an SGS, in 
which there is neghgible dispersion with respect to ky, the 
energy shift arising from the perturbation (or ky) should 
be neghgible. Assuming that k is small, the dominant 
summand is the case when s' — —s and m = 0. Thus, if 

{i^'-s k |cosa(a;) ay — sina(x) CTz] k) = , (11) 

the superlattice would be an SGS. Under this condition, 
the relative deviation of the energy dispersion relation. 



Eq. ([9]), from £'s(k) = s?ivo\kx\ (Eq. (3) in the paper) 
is O {ky/kxGx) and that of the wavefunction, Eq. (|10p . 
from Eq. ([8|) is 0{ky/Gx)- Similar quantities in a nor- 
mal graphene superlattice are O (ky/kl) and O {ky/kx), 
respectively. The eigenfunction ■0s.k(r) of the Hamilto- 
nian H in Eq. ^ is, to a good approximation, obtained 
by V's,k(r) = ■00 ski'") ^^'^ ^^^^^ eigenstate of ax 
because U commutes with a^- 



